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Abstract. Within the framework of ImQMD05, we study several isospin sensitive observables, 
• such as DR(n/p) ratios, isospin transport ratio (isospin diffusion), yield ratios for LCPs between 

| the projectile region and mid-rapidity region for the reaction systems Ni+Ni, Zn+Zn, Sn+Sn at 

£^ ■ low-intermediate energies. Our results show that those observables are sensitive to the density 

dependence of symmetry energy, and also depend on the cluster formation mechanism. By 
comparing these calculations to the data, the information of the symmetry energy and reaction 
^ ■ mechanism is obtained. 
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The nuclear symmetry energy plays an important role in the properties of nuclei and neutron 
stars [H El El [U [S] . To a good approximantion, it can be written as 

Esym = S( P )5 2 . (1) 

where 5 = {p n — Pp) / {Pn+ Pp), is the isospin asymmetry; p n , p p , are the neutron, proton densities, 
and S(p) describes the density dependence of the symmetry energy. Theoretical predictions for 
S(p) from microscopic nucleon-nucleon interactions show large uncertainties, especially in the 
region of suprasaturation density [SIE]. Constraining the density dependence of the symmetry 
energy has become one of the main goals in nuclear physics and has stimulated many theoretical 
and experimental studies El El [TO1 QH G21 H31 El H3 QS OH [El [13 EH ED E21 E3 El E3 ESI 
EH ESI [291 [30] . Heavy Ion Collisions (HIC) with asymmetric nuclei provide a unique opportunity 
for laboratory studies of the density dependence of the symmetry energy because a large range 
of densities can be momentarily achieved during HICs. In theoretical studies with transport 
models, the isospin ratio observables which are constructed from the isospin contents of emitted 
nucleons or fragments, such as Y(n)/Y(p) and DR(n/p) for emitted nucleons [201 EU [30], isospin 
transport ratios Ri constructed from the isospin asymmetry of projectile residues (or emitted 



resource) [TH [151 [261 [31] , an d -R^efcZ> constructed from the yields of LCP between the mid-rapidity 
and projectile region [32 J, have been prove to be primarily sensitive to the density dependence of 



the symmetry energy. By comparing the theoretical predictions to the experimental data, the 
sought-after constraints can be obtained. 

One frequently utilized transport models to describe the heavy ion collisions is the Boltzmann- 
Uehling-Uhlenbeck (BUU) equation, which provides an approximate Wigner transform of the 
one-body density matrix as its solution [33]. 
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The l.h.s. of this equation is the total differential of / with respect to the time assuming a 
potential U. Usually a Skyrme-parametrization of the real part of the G-matrix or Skyrme-like 
energy density functional are employed as the nucleonic potential which describe the influence 
of the different isospin asymmetric nuclear equation of state (asy-EOS). Stochastic extensions 
of these mean- field based approaches have been introduced (see [31] and references therein). 
For instance, in the so-called Stochastic Mean Field (SMF) [35] model fluctuations are injected 
in coordinate space by agitating the spacial density profile. The r.h.s. of Eq. ^fy contains a 
Boltzmann collision integral which describes the influence of binary hard-core collisions and are 
realized by the test particles. 

Another frequently utilized approaches, known as the Molecular Dynamics Model (QMD)that 
represent the individual nucleons as Gaussian "wave-packet" with mean values that move in 
according the Ehrenfest theorem; i.e. Hamilton's equations |36j. 
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Pi = -Vr,E(^) = -Vr,E( ff ) ( 4 ) 
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The expectation of the total Hamiltonian < H > is obtained from the real part of the G-matrix 
or the Skyrme energy density functional, and it describes the influence of the different asy-EOS. 
The collision part in the QMD models are handled as same as the way in BUU type models but 
it is for nucleons rather than the test particles. 

In this work, we choose to simulate nuclear collisions with the code ImQMD05 developed 
at the China Institute of Atomic Energy (CIAE), details of this code are described in Ref. 
[2H1 [301 EH [38] , for studying several isospin ratio observables, such as DR(n/p), isospin transport 
ratios E4, R™Jd ratios for the yields of LCP between the mid-rapidity and projectile region and 
their relations to the fragmentation mechanism. For brevity, we limit our discussion here to the 
parameterization of the symmetry energy used in our calculations, which is of the form 

s( P ) = \fpl / \^f 3 + ^r. (5) 

3 2m 2 po 2 po 

where m is the nucleon mass and the symmetry coefficient C s = 35.19MeU. Using this particular 
parameterization, the symmetry energy at subsaturation densities increases with decreasing 7$, 
while the opposite is true for supranormal densities. In general, the EoS is labeled as stiff-asy for 
7i > 1, and as soft-asy for 7$ < 1. Finally, we also give a brief discussion the recent comparisons 
between the ImQMD05 and SMF calculations for further understanding the theoretical issue in 
the describing the reaction mechanism at low-intermediate energy heavy ion collisions. 



2. isospin ratio from nucleon to fragments 

Isospin ratios which are constructed from the isospin contents of fragments, such as R(n/p) = 
Y(n)/Y{p) (or named as n/p ratio), DR(n/p) from neutron-rich and neutron-poor systems, 
isospin transport ratios R { and B$$ d = 2Y LC p(y° < 0.5)/Y LC p(0.5< y° <1.5), are sensitive to 
the density dependence of symmetry energy. In this section, we will check their sensitivities to 
the density dependence of symmetry energy and try to get the information of symmetry energy 
from them. 

2.1. n/p ratio and DR(n/p) ratio 

The neutron to proton ratio R n / P = Y{n)/Y{p) of pre-equilibrium emitted neutron over 
proton spectra was considered as a sensitive observable to the density dependence of symmetry 
energy [21], because it has a straightforward link to the symmetry energy. In order to reduce 
the sensitivity to uncertainties in the neutron detection efficiencies and sensitivity to relative 
uncertainties in energy calibrations of neutrons and protons, the double ratio 

DR(n/p) = R n/p (A)/R n/p {B) (6) 

had been measured by Famiano and compared with the transport model prediction [201 [2T] . 

We performed calculations of collisions at an impact parameter of b = 2fm at an incident 
energy of 50MeV per nucleon for two systems: A = 124 Sn + 124 Sn and B = 112 Sn + 112 Sn 
with ImQMD05 to study the DR(n/p) ratio for emitted nucleons [30] . The shaded regions in the 
left panel of FigQ]show the range, determined by uncertainties in the simulations, of predicted 
double ratios DR(n/p) = R n / p {A) / R n j p {B) of the nucleons emitted between 70° and 110° in 
the center of mass frame as a function of the center of mass nucleon energy, for 7, = 0.5 and 2.0. 
The double ratios DR{n/p) are higher for the EOS with the weaker symmetry energy density 
dependence ji = 0.5 than that for 7, = 2.0 because the nucleons mainly emit from the lower 
density region at intermediate energy HICs. Compare to the data on DR(n/p) for emitted 
nucleons (solid stars), the general trend of data DR{n/p) are qualitatively reproduced and the 
data seem to be closer to the calculation employing the EOS with 7, = 0.5. The right panel 
of figH] show the coalescence-invariant double ratio. The coalescence-invariant double ratios 
are constructed by including all neutrons and protons emitted at a given velocity, regardless of 
whether they are emitted free or within a cluster. The data are shown as open stars and the 
calculation results are shown as shaded region. Here, the measurement and simulation results 
illustrate that the fragments with Z > 2 mainly contribute to the low energy spectra and do 
not affect the high-energy DR{n/p) data very much. 

In order to constrain the range of 7$ from the DR(n/p) data that had been published, a 
series calculations for two systems, A = 124 Sn + 124 Sn and B = 112 Sn + 112 Sn, have been 
performed by varying 7$ = 0.35, 0.5, 0.75, 1.0 and 2.0[26]. Since the emitted nucleons are mainly 
from the subnormal densities at this energies, the n/p ratios of emitted nucleons are associated 
with the values of symmetry energy at subnormal density. Therefor, the DR{n/p) ratio should 
increase with decreasing 7$. However, in the limit of very small 7$ <C 0.35, the finite system 
completely disintegrates and the DR{n/p) ratio decrease and approaches the limit of reaction 
system, (N '/ 'Z)\^j '(N / 'Z)n2 = 1.2. As a consequence of these two competing effects, the double 
ratio values peak around the 7« = 0.7. Despite the large experiment uncertainties for higher 
energy data, those comparisons definitely rule out very soft (7, = 0.35) and very stiff (7, = 2.0) 
density dependence of symmetry energy. The x 2 analysis suggest that within a 2a uncertainty, 
parameters of 7$ fall in the range of 0.4 < 7^ < 1.05 for the C s = 35.2MeV . 

2.2. isospin transport ratio 

When the projectile and target nuclei come into contact, there can be exchange of nucleons 
between them. If the neutron to proton ratios of the projectile and target differ greatly, the 
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Figure 1. (Color online) (Left) DR(n/p) ratios for emitted free nucleons and (Right) coalescent- 
invariant DR{n/p) from the ImQMD simulations are plotted as shadow region. 



net nucleon flux can cause a diffusion of the asymmetry 5 reducing the difference between the 
asymmetries of two nuclei. This isospin diffusion process, which depends on the magnitude 
of the symmetry energy, affects the isospin asymmetry of the projectile and target residues in 
peripheral HICs. The isospin transport ratio Ri has been introduced [13] to quantify the isospin 
diffusion effects, 

2X — X aa — X bb 

Ki = V Y ' 

^-aa — ^-bb 

where X is an isospin observable and the subscripts a and b represent the neutron rich and 
neutron-poor nuclei. In this work, we use a and b to denote the projectile (first index) and 
target (second index) combination, where a = 124 Sn, and b = 112 Sn. We obtain the value of 
Ri by comparing three reaction systems, a + a, b + b and a + b (or b + a). Construction of 
the transport ratio minimize the influence of other effects besides isospin diffusion effects on 
the fragment yields, such as preequilibrium emission and secondary decay, by rescaling the 
observable X for the asymmetric a+b system by its values for the neutron-rich and neutron- 
deficient symmetric systems, which do not experience isospin diffusion. Based on Eq. ([7]), one 
expects Ri = ±1 in the absence of isospin diffusion and Ri ~ if isospin equilibrium is achieved. 
Eq. (|7|) also dictates that two different observables, X, will give the same results if they are 
linearly related. In one experiment, X was taken as the isoscaling parameter, a, obtained from 
the yield of the light particles near the projectile rapidity [12], to measure the isospin diffusion 
ability in heavy ion collisions. In transport models [14] [15], the isospin asymmetry 5 of the 
projectile residues (emitting source) has been used to compute Ri(5) because it is linearly related 
to the isoscaling parameters a[18l fl"2| 155]. 

We analyze the amount of isospin diffusion with ImQMD05 by constructing a tracer from 
the isospin asymmetry of all emitted nucleons (N) and fragments (frag), including the heavy 
residue if it exists, with velocity cut v^^ rag > 0.5v£ e ™ m (nearly identical results are obtained 
with higher velocity cut v^'^ ra9 > 0.7t^ e ™ m ). This represents the full projectile-like emitting 
source, and should be comparable to what has been measured in experiments. Fig. 2 shows the 
results of isospin transport ratios Ri(X = S^jrag) (upright triangles) as a function of the impact 
parameter for a soft symmetry case (7$ = 0.5, open symbols) and a stiff symmetry case (7, = 2.0, 
closed symbols). Ri obtained with soft-symmetry case is smaller than those obtained with stiff- 
symmetry potential case. This is consistent with the expectation that higher symmetry energy 
at subnormal density leads to larger isospin diffusion effects (smaller Ri values). 
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Figure 2. (Color online) Isospin transport ratios as a function of impact parameter with two 
tracers for a soft symmetry case (7, = 0.5, open symbols) and a stiff symmetry case (7, = 2.0, 
closed symbols). Upright triangle symbols are for the tracer defined by the isospin asymmetry 
of all fragments and unbound nucleons with velocity cut {v^^ ra9 > 0.5t>g e ™ m ), X = SNjrag- 
Circles are for the tracer defined by the heaviest fragment with Z max > 20 in projectile region, 

X = $Z max >20- 



Ri depends weakly on impact parameter over a range extending from central (b = 3fm) to 
mid peripheral(6 = 8fm) collisions. Interestingly, the isospin equilibrium and global thermal 
equilibrium are not reached even for central collisions. Our results show, that neither the 
effective interaction is sufficiently strong nor the collisions are sufficiently frequent (most of them 
are Pauli suppressed) to mix the projectile and target nucleons completely. These two effects 
prevent the combined system from attaining isospin equilibrium even in central collisions. With 
impact parameter increasing for b > 5fm, the overlap region and thus the number of nucleons 
transferred from projectile and target decreases, causing the Ri values to increase. 

In peripheral collisions, most often, a large residue remains. If it decouples from the full 
emitting source before it equilibrates, it may experience a different amount of diffusion than the 
full emitting source examined by X = 5 n, frag- To examine this, we constructed a tracer using 
the isospin asymmetry of the heaviest fragments with charge Z max > 20 in the projectile region. 
This tracer is mainly relevant to peripheral collisions as the central collisions are dominated 
by multifragmentation and very few large projectile fragments survive. The dependence of 
Ri{X = Sz max >2o) for impact parameter b > 5/m is shown as open and closed circles in Fig. 
2. The isospin transport ratios constructed from the different isospin tracers have different 
values especially in the case of ji = 0.5. Stronger isospin equilibration (smaller Ri values) 
is observed in the isospin transport ratios Ri(X = dj^jrag) constructed from nucleons and 
fragments than Ri(X = 5z max >2o) constructed from the heaviest fragments with Z max > 20. 
Since isospin diffusion mainly occurs through the low-density neck region, and the system 
breaks up before isospin equilibrium, the asymmetry of the projectile and target residues do 
not achieve equilibrium and, larger Ri(X = dz max >2o) values result. In contrast, there is more 
mixing of nucleons from the target and projectile in the neck region due to the isospin diffusion. 
Consequently, rupture of the neutron-rich neck is predicted to result in the production of neutron- 
rich fragments at mid rapidity. 

Since fragments are formed at all rapidities, we can examine the rapidity dependence of Ri 
to obtain more information about the reaction dynamics. Fig. 3 shows Ri as a function of 
the scaled rapidity y/ubeam- The symbols in the leftmost panel are experimental data obtained 
in Ref. [18] for three centrality gates. This transport ratio was generated using the isospin 



tracer X = ln(Y( 7 Li)/Y( 7 Be)) where Y( 7 Li)/Y( 7 Be) is the yield ratio of the mirror nuclei, 
7 Li and 7 Be [18J. As expected the values of Ri obtained from peripheral collisions (solid stars) 
are larger than those obtained in central collisions (open stars). For comparison, the ImQMD05 
calculations of Ri(X = 5 at, frag) are plotted as lines in the middle and right panels for a range 
of impact parameters. The middle panel contains the results from the soft symmetry potential 
(li = 0.5) while the right panel shows the results from the stiff symmetry potential (7, = 2.0). 
The impact parameter trends and magnitude of the data are more similar to the results of the 
calculations from soft symmetry potentials (7$ = 0.5) for peripheral collisions. 

We have performed the \ 2 analysis for both observables, Ri and Ri(y), for constraining the 
density dependence of symmetry energy. Using the same 2a criterion, the analysis brackets 
the regions 0.45 < 7$ < 0.95 is obtained. It is consistent with previous analysis on DR(n/p). 
However, the experimental trend of Ri gated on the most central collisions (open stars) is 
not reproduced by the calculations. The experimental data indicate more equilibration for 
central collisions near mid rapidity while the transport model indicates more transparency. The 
equilibration in the E /A = 50MeV data may be the result of the impact parameter determination 
from charged particle multiplicity wherein the most central collisions are assumed to be the ones 
with highest charge particle multiplicity. For the most central events, a gate on the highest 
multiplicity, may select events in which more nucleon-nucleon collisions occur rather than a 
strict selection on the most central impact parameters. 
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Figure 3. (Color online) (Left panel) Experimental Ri as a function of rapidity for three 
centrality gates [16]. (Middle panel) The calculated results of Ri(X = 5^j rag ) as a function of 
rapidity for 6 = 2,4, 6, 8fm for 7$ = 0.5 and (Right panel) 7$ = 2.0. 



2. 3. R™i l Jd ratios for light charged particles 

The yield ratios, R^fJ^, for LCPs between the projectile region and mid-rapidity region, are 
defined as 



ymd 2-Yield(0.0<Y r <0.5) 



"-yield- Yield(0.5 < Y r < 1.5) ' { ) 

where Y r = is the reduced rapidity. It reflect the isospin migration ability and have been 

proj 

measured by Kohley.et.al. for 70 Zn + 70 Zn, 64 Zn + 64 Zn, 64 Ni + 64 Ni at the beam energy of 35 
MeV/nucleon for middle peripheral collisions [32] . The data show a clear preference for emission 
around the mid-rapidity region for more neutron-rich LCPs resulting from the isospin migration 
mechanism through the neck region between the projectile and target [4*0|I29|. Theoretical study 
by SMF model[H] demonstrates that R™£j[d ls sensitive to the density dependence of symmetry 



energy. The experimental trends are reproduced by SMF model. However, there are largest 
discrepancies on the reduced rapidity distribution for the yields of proton and 3 He, and also 
on their values of -R^fd as mentioned in Reference |32j. The discrepancies were thought to be 
related to the statistical decay of QP at later stages of the reaction [32J . In the points of reaction 
dynamics, the different fragmentation mechanism in the transport models simulations also lead 
to different behaviors of the rapidity distribution for LCPs besides the effects from secondary 
decay. Thus, it is instructive to study the rapidity distribution of LCP with ImQMD05. 

In FigHKa), we present the multiplicity distribution for fragments with Z > 3 for 70 Zn + 70 Zn 
at Eb eam =35 MeV/nucleon and impact parameter b=4fm and 73 = 2.0. We find that half of 
events belongs to multi-fragmentation process which are defined by multiplicity for fragments 
with charge Z > 3, i.e.,M(Z > 3) > 3. The rest are the binary (M{Z > 3) = 2) and ternary 
{M(Z > 3) = 3) fragmentation events. It suggests that the binary, ternary fragmentation and 
multi-fragmentation coexist around 35MeV/nucleon. In FigJH (b) and (c), we plot the reduced 
rapidity distribution for the yields of 3 He and 6 He obtained with three kinds of fragmentation 
process, binary (square symbols), ternary (circle symbols) and multi- fragmentation (triangle 
symbols) which are selected by M(Z > 3) = 2, 3 and > 3. The yields of s He and G He in FigJ5] 
are normalized to per event. It is clear that the binary events produce more 3 He and e He 
at mid-rapidity relative to that produce in multi-fragmentation events. For 7$ = 2.0 case, the 
yield of 3 He at Y r = obtained with binary fragmentation events is 35% larger than that with 
multi-fragmentation events. For neutron-rich LCP, for example, the yield of e He at Y r = 
obtained in binary fragmentation events is 70% larger than that in multi-fragmentation events 
due to isospin migration. 
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Figure 4. (Color online) (a) The multiplicity distribution for fragments with Z > 3 {M(Z > 3)). 
(b) is the reduced rapidity (Y r ) distribution for the yield of 3 He with binary (square symbols), 
ternary (circle symbols) and multi-fragmentation (triangle symbols) process, (c) is for G He. All 
of those results are for 70 Zn + 70 Zn at E=35 MeV/u for b=4fm and ji = 2.0. 



Figj5] shows the calculated results for the rapidity distribution of light charged particles p, 
d, t, 3 He, A He and 6 He for m Ni + 64 Ni at b=4fm with 100,000 events. The distributions are 
normalized with the yield at Y r = for comparing with data in Ref [32] to understand the 
fragmentation mechanism, because the fragmentation mechanism mainly determine the shape 
of the rapidity distribution for LCP. The open circles are for 7$ = 0.5 and solid symbols are 
for 7j = 2.0. The data are taken from the Ref [32] and plotted as stars. Both our calculations 
and data show that the width of distribution decreases with the mass of LCPs increasing. For 
the rapidity distributions of 3 H and 3 He, the width of distribution for 3 H is smaller than that 



for 3 He due to the isospin migration. By comparing the simulated results to the data, the 
ImQMD05 calculations with stiffer symmetry energy can well reproduce the data at forward 
rapidity region (Y r > 0) for all p, d, t, s He, 4 He and e He. For the backward rapidity region 
(Y r < 0), there are obvious differences between the results from ImQMD05 calculations and the 
data because the efficiency for detection of LCPs at the backward [32] are not included in this 
ImQMD05 calculations. 
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Figure 5. (Color online) Reduced rapidity (Y r ) distributions for p, d, t, 3 He, 4 He and e He 
fragments from the 35 MeV/nucleon 64 Ni + 64 Ni reaction for impact parameter b=4fm. The 
experimental data are shown as the stars. The ImQMD05 calculations for 7, = 0.5 are shown 
as the open circles and the solid circles are for 7^ = 2.0. Each distribution is normalized with 
the yield at Y r = 0. 



In order to constrain the symmetry energy by the rapidity distribution of LCPs, we further 
calculate K^fid m a ser i es of 7« = 0.5, 0.75, 1.0, 2.0. In Figj6j we present the results of K^fid as a 
function of the AZ of emitted particles for three reaction systems 64 Zn + 64 Zn, 64 Ni + 64 Ni and 
70 Zn + 70 Zn at b=4fm. The open symbols are the results for 7, = 0.5, 0.75, 1.0 and 2.0. The solid 
stars are the data from [32]. Since the isospin migration occurs in the neutron-rich neck region, 
the -R^gfd shows an increasing trend with the values of isospin asymmetry of LCPs increasing for 
the same element. The more neutron-rich the LCPs is, the larger the -R^fd Furthermore, the 
calculated results show the values of -R^fd for neutron-rich isotopes are sensitive to the density 
dependence of symmetry energy. The calculations with stiffer symmetry energy predict larger 
values of due to the stronger isospin migration effects. This conclusion is as same as the 

results obtained with SMF model [32J. As shown in Figj6l the ImQMD05 calculations with stiffer 
symmetry energy (7, > 0.75) can reasonably reproduce the data of -R^fd as a function of AZ for 
64 Zn + 64 Zn. But our calculations underestimate the R^fJ^ values of neutron-rich light charged 
particles, such as 6 He, for the neutron-rich reaction systems 64 Ni + 64 Ni and 70 Zn + 70 Zn. It 
could come from the lacking of fine structure effects of neutron-rich elements (such as neutron- 
skin, stability of lighter neutron-rich elements), and the impact parameter smearing effects in 

ymid 
yield 
imid 
yield 

l yield 

reaction system. 



the transport model simulations. Even though our calculations can reproduce the -R^w data for 
Zn + 64 Zn, the definitely constraints on the symmetry energy with the data of R^ ld can not 
be obtained before we fix the problems on the theoretical predictions of -R^llfd f° r neutron-rich 
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Figure 6. (Color online)i?^^ rf values as a function of the charge times mass (ZA) for p (ZA=1), 
d (ZA=2), t (ZA=3), 3 i?e (ZA=6), 4 #e (ZA=8), 6 He (ZA=12). The open symbols are the 
results obtained with ImQMD05 for 7j = 0.5, 0.75, 1.0 and 2.0. The solid stars are the data from 
[32]. 



2-4- remarks on the comparisons between ImQMD05 and SMF 

Up to date, there are several constraints on the symmetry energy with isospin sensitive 
observables, DR(n/p) raitos, isospin transport ratios Ri and R^jfj ratios, by adopting different 
type of transport models, such as QMD type and BUU tvpepUl H3J HI] . There are overlap 
between the results of symmetry energy from different approaches, but they are different in 
detail. Thus, further understanding the issues in the QMD type and BUU type would be crucial 
in theoretical studies for improving the constraints on the symmetry energy. 

At the code level, both BUU and QMD models propagate particles classically under the 
influence of a mean field potential, which is calculated self-consistently the positions and 
momenta of the particles, and allow scattering by nucleon-nucleon collisions due to the residual 
interaction. The Pauli principle in both approaches is enforced by application of Pauli blocking 
factors. These similarities in implementation have lead to similarities in predictions for many 
collision observables [45j . 

There are also significant differences in these approaches. In the BUU equations, each nucleon 
is represented by 200-1000 test particels that generate the mean field and suffer the collisions. 
In QMD, there is one test particle per nucleon. A-body correlations and cluster formation are 
not native to the original BUU approach; which is supposed to provide the Wigner transform of 
the one body density matrix. On the other hand, many-body correlations and fluctuations can 
arise from the A-body dynamics of QMD approach. Such A-body correlations are suppressed 
in BUU approach, but correlations can arise in both approaches from the amplification of mean 
field instabilities in spinodal region [M]. Collision algorithms in the QMD approach modify the 
momenta of individual nucleons, while in BUU approach, only the momenta of test particles 
are modified. Depending on the details of the in-medium cross sections that are implemented, 
the blocking of collisions can also be more restrictive for QMD than for BUU, leading to fewer 
collisions and therefore a greater transparency. 

Fragments can be formed in QMD approaches due to the A-body correlations and these 
correlations are mapped onto the asymptotic final fragments by a minimum spanning tree 
algorithm. Serval different methods have been developed to allow BUU codes to calculate cluster 
production. In the Stochastic Mean Field (SMF) approach (one of the BUU type model), the 
time evolution of the one-body phase-space distribution / is governed by the nuclear mean- 



field, two-body scattering, and a fluctuating (stochastic) term which causes the fragmentation 
|35l W\\ I44j , Since there are typically more than 100 test particles per nucleon, collision 
induced fluctuations are smaller in BUU than in QMD possibly suppressing the fragment 
formation rates. 

As an example, we present the results of average charge number for Z > 3 as a function of 
rapidity obtained with ImQMD05 and SMF models for 124 Sri + 124 Sn at E beam = 50AMeV 
and b=6, 8fm in Fig0 The solid lines are the results from ImQMD05, and dashed lines are 
from SMF calculations [42 j . Squares are for b=6fm, and circles are for b=8fm. The results 
obtained with QMD simulations show that the average charge number of Z > 3 increase with 
rapidity increasing before y/y^ m ~ 1 in the forward region in the QMD simulations, and 
two peaks appear around the projectile and target rapidity region, respectively. In the SMF 
simulations, the peak appears at mid-rapidity besides two peaks around projectile and target 
rapidity. It clearly shown that the strict Pauli blocking in the QMD models simulations leads to a 
greater transparency than that in the SMF type simulations. The larger fluctuation in the QMD 
models lead to more fragments and light particles emitted than that from SMF predictions. As 
results, the fragments distributed over the whole rapidity region and heavier fragment is, larger 
velocity is. In SMF calculations, the fragmentation process is less effective due to the reduced 
amplitude of fluctuations and many-body correlations. This enhances the appearance of binary 
and ternary processes in semi-peripheral heavy ion reactions, according to the prominent role of 
the mean-field dynamics. As a consequence, the intermediate mass fragments tend to distribute 
at mid-rapidity. It also lead that the the average charge number of Z > 3 have a peak around 
mid-rapidity and narrower rapidity distribution of lighter clusters than the results from QMD 
models. 




Figure 7. (Color online) The average charge number for Z > 3 as a function of rapidity for 
124 Sre _|_!24 g n a £ D= g 5 8f m w ith 73 = 2.0. The solid lines are the results from ImQMD05, and 
dashed lines are from SMF. 



3. Summary 

In summary, we have investigated the influences of the density dependence of the symmetry 
energy on several different isospin ratio observables, such as DR(n/p) ratio, isospin transport 
ratios Ri, the rapidity dependence of isospin transport ratio Ri(y) and -R^fd r &itos (the yield 
ratios of LCP between the mid-rapidity and projectile region) with ImQMD05. The study 
shows that these isospin ratio observables are sensitive to the density dependence of symmetry 
energy. This conclusion is similar to conclusions reached using BUU approaches in the range 
of symmetry energies studied here. By comparing the calculated results to data, the very soft 
(7i = 0.35) and very stiff symmetry energy (7$ = 2.0) are ruled out. 



Cluster formation is important for intermediate energy heavy ion collisions, and it modifies 
the spectral double ratios at E cm _ < AOMeV. We also tested different tracers by constructing 
corresponding isospin transport ratios for them using different symmetry energies. For weakly 
density dependent symmetry energies (small 7$) with large symmetry energies at sub-saturation 
densities, the values of Ri for the residue tracer X = 5z max >20 are larger than those extracted 
from the entire emitting source, i.e., X = 5 n, frag- The difference between these two tracers can 
be examined experimentally as a new probe of the symmetry energy and reaction mechanism. 

By studying reaction systems 64 Zn + 64 Zn, 64 Ni + 64 Ni and 70 Zn + 70 Zn at the beam energy 
of 35 MeV per nucleon and b=4fm within the framework of ImQMD05, we find that half of 
events belongs to the multi-fragmentation mechanism, and half of them is of binary and ternary 
fragmentation events. The binary and ternary events produce more light charged particles at 
middle rapidity, and the multi-fragmentation events broaden the reduced rapidity distribution 
for the yields of LCPs. Both the data and our calculations illustrate that the reaction systems 
seems more transparency and more fragments, light particles emitted. As results, the data of the 
reduced rapidity distribution for the yields of LCPs and R^fJ^ as a function of AZ for m Ni+ m Ni 
can be well reproduced by the ImQMD05 calculations. For neutron rich reaction systems 
64 Ni + 64 Ni and 70 Zn + 70 Zn, our calculations underestimate the -R^fd values of neutron-rich 
light charged particles, such as Q He, it could be cause by the lacking of fine structure effects for 
lighter elements in the transport models, and the impact parameter smearing effects. 
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